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Abstract. This article extends the non-extensive entropy of Tsallis and uses this 
entropy to model an energy producing system in an absorbing heat bath. This 
modified non-extensive entropy is superficially identical to the one proposed by 
Tsallis, but also incorporates a "hidden" parameter that provides greater flexi- 
bility for modeling energy constraints. This modified non-extensive entropy also 
leads to a more generalized family of energy transformation functions and also 
exhibits the structural scale invariance properties described in a previous article. 
This energy transformation also provides a more natural mechanism by which 
arbitrary power-law distributions can be stated in exponential form. 
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1 Introduction 

In a recent article, Fleischer 1 1 1 described a number of scale invariant and symme- 
try properties in non-extensive systems — systems based on a non-extensive form 
of entropy developed by Tsallis \2\. These scale invariant properties show that in 
aggregating energy states of a large ensemble, many of the mathematical relationships 
associated with a given state also apply to the aggregations and use identical mathemat- 
ical forms. Investigation of this scale invariance property revealed an energy transfor- 
mation function. The frequency of its appearance in these scale invariant mathematical 
forms suggests that energy transformations are an important element in modeling non- 
extensive systems. This energy landscape transformation also provides an additional 
perspective on non-extensivity. Tsallis' entropic form essentially shifts the probability 
weights in the system and thereby creating a power law distribution. This same effect 
can be captured by the energy landscape transformation functions. Thus, the energy 
landscape transformation function apparently lies at the heart of the scale invariance 
and symmetry properties present in non-extensive systems and seems to account for 
many of the useful applications of the Tsallis entropy especially in modeling dissipative 
systems where energy producing systems in a heat bath lose energy to its environment. 

In this article, new forms of scale invariance and energy transformations are illumi- 
nated by generalizing (extending) Tsallis' entropy formulation to explicitly account for 
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and generalize this energy transformation property. In this vein, Tsallis' entropy and all 
of its associated scale invariance and symmetry properties become special cases in a 
larger array of scale invariance and symmetry. These new properties show a number of 
features involving recursion (aggregations of aggregations) and ways of characterizing 
power-law distributions using exponential forms where the exponent is the transformed 
energy value. In this sense, this article expands on the results in Q. 

This article is organized as follows: Section[2]provides some background on Tsal- 
lis entropy and the scale invariant properties described earlier. Section|3]describes the 
new form of entropy and related system constraints. These two components are then 
used to define a new stationary probability. Section|4]briefly describes the scale invari- 
ance properties based on the forms of the new entropy and stationary probability and 
the associated/am/Zy of energy transformation functions. Section|5]explores the related 
symmetry and power law relationships. Finally, Section |6] offers some discussion on 
the implications of this modified Tsallis entropy, energy transformations and practical 
applications in the field of complex system simulation. Section ^provides some con- 
cluding remarks. 



2 Background 



Tsallis 1 3 2 1 developed a new entropy expression that forms the basis of a non-extensive 
form of thermodynamics: 
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S q = : '- (1) 

q-1 

where A: is a constant and S q is the entropy parameterized by the entropic parameter q. 
In classical statistical mechanics, entropy falls into a class of variables that are referred 
to as extensive because they scale with the size of the system. Intensive variables, such 
as temperature, do not scale with the size of the system. 1 Tsallis' form of entropy is 
non-extensive because the entropy of the union of two independent systems is not equal 
to the sum of the entropies of each system. That is, for independent systems A and B, 

S q (A + B)= S q {A) + S q (B) + ^-Q)S q (A)S q (B) _ ^ 

Tsallis uses this entropy to calculate a stationary probability pi for a canonical system 
where S q is maximized. Tsallis shows that pi is distributed by a power law 

^q ^ q 



1 Combine two vessels of gas each with the same volume and pressure into another vessel of 
twice the volume and the pressure and temperature of the combined gas will be the same as 
before. Energy and entropy, however, are examples of extensive variables in that combining 
several sources of either energy or entropy and you increase the total energy or entropy. 
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where a = (q—l)/t and based on system constraints (but see 1 1141 for a discussion on 
constraints) 

w 

X> = 1 (4) 

i=l 
W 

y^Pi/i — U, a constant. (5) 

i=l 

This distribution is different than the exponential law in the classic Boltzmann-Gibbs 
distribution, symbolized here by tt,. Tsallis points out that for q — > 1 the extensivity 
properties of classical statistical mechanics emerge. Thus, e.g., lim g ^i S q = Sbg = 
—k Y^,i lriTTi, an d lim g _,x = ""i(i)> nence me Tsallis entropy is a generalization 
of the Boltzmann-Gibbs entropy 12131 . 

Tsallis 1 5 1 also notes that the form of is the simplest form that satisfies certain 
non-additivity assumptions and uses only one parameter, namely the entropic parameter 
q . In the sections that follow, we both retain this inherent simplicity, but at the same 
time incorporate a new, "hidden" parameter m that provides additional flexibility for 
modeling dissipative systems. 



2.1 Scale Invariance and Symmetry in Extensive and Non-Extensive Systems 



In 1 1 6 1 a number of scale invariant properties in both classical and non-extensive sys- 
tems are described and are based on aggregations of energy states. Briefly, this type of 
structural scale invariance is best illustrated by the following: for any aggregated set of 
energy levels A = {ix, 12, ■ ■ ■ , in} where index ik refers to some particular energy level 
(note A can also simply consist of a single energy energy level i), 

d -^ = ^[fA(t)-(f)(t)] (6) 

where 

n A (t) =$>(*) and f A (t) = ^A*fh 

with the latter equivalent to the conditional expectation of a random energy value / 
given that the current state i is in set A. A similar form of scale invariance exists for 
second moments. Thus, 

df A (t) _ 

dt t 2 K) 

where o\ (t) is the conditional variance of energy (objective function) values at tem- 
perature t given that the current state being in set A. See |6, p. 232-33] for details and a 
more formal treatment. 

Fleischer 1 1 1 demonstrated a similar result for non-extensive systems where the 
only difference between the scale invariance properties of the classic case and that of 
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the non-extensive case was that the latter involved an energy transformation function. 
Thus, for any aggregated set of energy levels A, 



dp A (t) PA(t) 



dt 



t 2 



where 



p A {t)=J2Pi{t) and h{t) = 



fA(t) ~ (/>(*) 



where the transformed energy value /$ is defined by 

fi 



fi(fi,q,t) = 



PA® 



(8) 



(9) 



(10) 



where a = (q — l)/t. It will henceforth be notationally convenient to simply refer to 
transformed energy value using fa (t) , fa or f a as the case may be where it is understood 
to depend on the parameters i, fi, q, t. 

Fleischer 1 1 1 also showed that the scale invariance in second moments incorporates 
this energy transformation, hence for all aggregations A 



djA 
dt 



(f 2 )A-{f}A + (t-m 2 ) A 

t 2 

*k , (9-1)(/ 2 )a 
t 2 t 2 



(11) 



2.2 Discussion 

The appearance of the transformed energy value / in these scale invariant forms sug- 
gests that it is an important component in non-extensive statistical mechanics. It is 
therefore reasonable to infer that for an energy absorbing heat bath, the "equilibrium" 
condition associated with the internal energy constraint in l|5} implicitly involves this 
transformed energy value. One viewpoint suggests that in an energy producing system 
in thermal equilibrium with an energy absorbing heat bath there must be some way of 
characterizing the rate at which the energy is produced and absorbed. To generalize this 
notion further, a heat bath that absorbs energy at a high rate, must be balanced by a 
higher rate (in some sense) of energy production ;/ one seeks to model some equilib- 
rium condition. We can then perhaps model different rates of energy production and 
absorption by making the absorption rate proportional to some function of the energy 
level in a canonical system. In this way, the notion of non-extensivity can be further 
expanded and modeled in a fashion that encompasses the information/entropy loss, as 
indicated in for q > 1, and an equilibrium condition that encompasses an energy 
production/absorption component. 

To capture these notions in a more flexible way, a modified Tsallis entropy S u is de- 
fined and in a manner that provides additional flexibility for modeling canonical systems 
in energy absorbing heat baths yet retains the inherent simplicity of the non-extensive 
entropy. The increased flexibility is based on using a parameter u that involves both the 
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entropic parameter q from Tsallis' entropy and an energy parameter m. Note that in 
Tsallis' approach, the entropic parameter q is typically associated with the state proba- 
bility pi (see e.g., Q and Tsallis notes the effects of the exponent q in the Type 2 
constraint as they "[privilegiate] the rare and the frequent events" depending on whether 
q < 1 or q > 1, respectively Q p.535]. But this notion of shifting the probability weight 
of different energy values is, in some sense, equivalent to transforming the energy val- 
ues themselves. Thus, the entropic parameter q plays a central role, in capturing this 
shift of probability by transforming the energy landscape as indicated in the derivation 
of the energy transformation function in i25\ . Thus, q > 1 implies that the exponent 
of fi in is also 1 . It would however be useful to permit some additional freedom 
in choosing the exponent of fi to capture this notion of an energy producing system in 
equilibrium with an absorbing heat bath. 

Taking these considerations into account, the parameter u which involves both q 
and m should be associated with the state probability pi, and so will also affect the 
energy transformation function. Because we also want to fashion situations where the 
exponent of fi in is not always equal to 1, we require that the exponent of fi in (0 
be different from that of pi. Keeping these qualifications in mind, the following section 
describes a modified non-extensvie entropy S u where a number of new scale invariant 
properties emerge. 



3 The Modified Non-Extensive Entropy S u 

In this section, Tsallis' entropy is generalized to increase the flexibility in using energy 
transformation functions. In this formulation, the Tsallis entropic parameter q in ([0 is 
replaced by an entropic/energy parameter u (u always follows q) that simultaneously 
accounts for both the parameter q and an energy transformation parameter m > 1 where 

u=(g-l)(m-l) + l (12) 

and the modified Tsallis entropy is defined by 



u - 1 



(13) 



Note that this form is identical to Q, hence retains all of its inherent benefits. It is easy 
to see that 

Km S u = lim S u = lim S q = Si = S Ba (defined earlier) (14) 

u — >1 q — >1 q — >1 

and certain other relationships and properties are easily extended and generalized. For 
example, (0 becomes 

S U (A + B) = S U (A) + S U (B) + (l-^(A)S u (B) 

= S U (A) + S U (B) + ^-^-^(A)S U (B) 

k 
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since u — 1 = (q — l)(m — 1). Thus, for any given value of g > 1 where some entropy 
loss occurs, the magnitude of this loss can also be modeled using the parameter m 
which, as explained below, is associated with the energy levels. Note that if q = 1, the 
value of m becomes irrelevant. The next section explores the implications of this simple 
modification. 



3.1 The Energy Loss Rate 

The parameter m is useful for modeling an energy producing system in thermal equilib- 
rium with an absorbing heat bath. Such a system dissipates its energy to its surroundings 
while maintaining an average value of its internal energy. A canonical ensemble of such 
a system can be modeled therefore by a power of the energy function / as in 

w 

^Pifi"^ 1 = Ua. constant (16) 

i=l 

where the parameter m serves to capture the notion of an energy loss rate or energy 
dissipation rate (or perhaps, the energy absorption rate). The exponent u of the station- 
ary probability pi thus serves the same purpose as in Tsallis' works (see B4I7I Notice 
that (rn = 2) => (u = q), and the resulting system is equivalent to those based strictly 
on S q . Notice also that for m = 2, this constraint has the same form as that described 
by Tsallis' Type 2 constraint in (|5j and also in Tsallis' Type 3 constraint involving the 
so-called "escort probabilities" (see |7 1). 



3.2 The Stationary Probability 

Tsallis 1 2 7 1 illustrates how maximizing the entropy S q given the constraints in and 
leads to a stationary probability and results in the well-known stationary probability 
in Q. To obtain the stationary probability pi in light of the modified Tsallis entropy S u 
subject to the normalization constraint in (0} and the internal energy constraint in d!6i 
we use a similar approach as in |2 4|. The general Lagrangian function is 



L = k\ 
Therefore, 



TZiPi 



w 



5>(*)-i )-0 nr^/r- 1 -^ 



dpi 



U Pi 



„«-l 



u— 1 r m — 1 



1 



U 



+ a-pu P r 7, 



Setting Jl 81 to zero and rearranging, we obtain 

[i + Piu-i)!™- 1 } 



1 



and hence 



a{u — 1) 



[l + /3(«-l)/ ( 



m— 11 u-i 



(17) 



(18) 



[1 + (3(u - 1)/, 



m— 11 l- 



(19) 
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where f3 is often symbolized by the inverse temperature 1/t and 

Z u = ^[l+/3(u-l)^- 1 ]^ 

is the corresponding normalization constant. Note that dl9l l has the same form as Q 
except for the exponent of /,;. Noting that u—1 = (q — 1) (m — 1) we obtain the general 
form 

_ [l + Z^-lXm-l)//"- 1 ] 



(l-,)(m-l) 



Letting a = — 1) = (g — l)/t for notational convenience, then 

[1 + o(m - l)/™" 11 



(20) 



Notice that for the case m = 2, (I20> is equivalent to l|3}. 



4 SA Scale Invariance Based on S,, 



To demonstrate scale invariance based on aggregated states in the non-extensive case 
involving S u , we proceed in a similar fashion as in [ 1 1 while making the necessary 
adjustments to account for the definition of u. The following theorem states a scale 
invariant mathematical structure for systems based on S u . 

Theorem 1. Let A = . . . , i n } be any aggregation of energy levels i where the 

energy level associated with i is denoted as /j. Using the definition ofpi in i !2Qt . define 
the stationary probability of set A by pa — SieA Pi an d energy value associated with 
aggregated sets A by 

Tlii£A Pi ft 



f ' 

J A 



PA 



where 



P n 
J i 



fl 



l + o(m-l)/; 
Define the conditional variance of f A TL ~ 1 by 



I i>2m — 2\ I 2m— 1\2 

= \J )A - (J ) A - 



Then 



and 



dpA 
dt 



dt 



Pa 
t 2 



t 2 



£ m—1 I £i 
J A ~ \J 



- 1 ) 



{u-l){f 2m - 2 ) J 



8 



Proof: 

It is convenient to first consider the relevant quantities associated with individual states 
(energy levels) i, which in this case refers to dpi(t)/dt and d(f i m ~ 1 )/dt where these 
quantities are defined below. We proceed in similar fashion as in 1 1161 taking into ac- 
count the definition of u, S u and the exponent of /j. 

For notational convenience and simplicity, let Ni(t) be the numerator in J20I . Thus, 



Ni(t) = 



1 



u-l 
t 



n 



m — l 



(21) 



(hereinafter we will drop the (t) from Ni (i) to further simplify the expressions) and 
taking the derivative of i20\ with respect to temperature t, 



dPi(t) _ Z u §- t N l -N l §- t Z u 



dt 



In this case, 



(22) 



ON, _ d \, (u-l 
~dt ~ dt 

with Z u =J2j Nj, and Pi(t) — Ni/Z u and hence 



t 2 



(23) 



dt 



3 3 

Substituting this and d23i into (I22> yields 



dt 



t 2 



at 



(E.7 N 3? 
ju — 1 £ m—l 



J 3 



r 2 



Pi® 
t 2 



atu-1 pm-1 _ ^3 3 J 3 
■i J i 



m — l 



To further simplify the notation, define the transformed energy value 

p m — l _ t\ju — 1 r m — 1 
J i — ^ i J i 



V 



n 



l + a(m-l)// 



Mi. 



where a = (q — l)/t and, 



(24) 



(25) 



Mi = 1 + 0(70-1)/; 



(26) 
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(the reason for this latter definition will become clear later on) without the clutter of the 
arguments. Substituting d25l into d24b and further simplifying yields the two equivalent 
forms: 



dpi 
dt 



Pi 
t 2 
Pi 
t 2 



Mi 



if™' 1 ) 

fm-l 

M 



Now, taking the derivative of the probability of the aggregated set, 

dp a d sr^ dp* 



dt dt 



dt 



(27) 
(28) 

(29) 



ieA ieA 

and substituting &27\ into ( I29l l (keeping in mind the dependence on t) yields 

dp 



"PA _ Pi_ [fm-l _ ilm-X\ 
dt ^ t 2 l Jl U 



ieA 



E 

ieA 



PA 



Pif? 



t 2 

Pif, 



E 



< 2 £ : PA 



Piif" 1 - 1 ) 
t 2 

PAif" 1 - 1 ) 



i£A 

TO— 1 „ /i'm— 1 



t 2 



PA 
t 2 



- if™- 1 ) 



(30) 



where for aggregated states, J30I has a similar mathematical structure as J27I28> . hence 
exhibits a scale invariance property the foundation of which is based on the energy 
transformation function f™~ . 

Scale invariance for second moments is indicated in the following where again the 
energy transformation of /,; is used. Keeping in mind the dependence of / and pi on the 
temperature t, consider 



d(f m ^) _ df : 



m— 1 

a 



dt 



dt 



d_ 

dt 



Uen 



ien 



E 

ien . 



dPi fra-l 

dt h 



dfl 



dt 



(31) 



Substituting 1271 into the first part of (13 1 1 and simplifying yields 



dt 



ien 



ien 



t 2 



ien 



dfl 



dt 



(32) 
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Noting the form of the first two terms on the right-hand-side in i32\ and the fact that in 
the third term 

dfi" 1 _ }2m-2 f u ~ M nov 
dt h \ t 2 ) Kii) 

and substituting into fl32i and adding the symbol J? to denote expectations over the 
entire state space yields 

d(/ m -% _ (/ 2r - 2 )r2 - (t^fn + (" - l)(/ 2m ' 2 )r2 
St t 2 

1 (34) 



?h , («-i)(/ 2m - 2 ), 



i 2 i 2 

where a 2 n represents the variance of the values //" over the entire energy landscape 
(at temperature t). 

Eq. d34l provides the basis for another form of scale invariance. Thus, after going 
through similar steps as in (13 1 1 through (134ft we get 



dt t 2 PA t 2 \ t 2 J PA 



(35) 



Noting that the first and third terms indicate conditional expectations conditioned on 
the current state being in set A, then d35ft can be re-written in the convenient notation 



dt t 2 

- IT + ^ ( 36 ) 

where d36i is clearly analogous to (134i and so exhibits a form of scale invariance. Note 
that the terms involving (it — 1) also scale with the aggregated set A in the non-extensive 
case. See 111 |6| for similar results in the case of S q and Si, respectively. ■ 

An interesting aspect of these relationships can be succinctly described using the 
following commutative-like diagram in Figure^ A similar diagram can be produced to 
depict an analogous relationship for second-moments. The next section further explores 
some of these symmetry relationships. 



5 Symmetry Relationships 

Fleischer [ 1 1 explored a number of symmetry relationships in addition to ones Tsallis 
has indicated in his early works based on the effects of different constraints used in 
modeling a canonical system (see also [7|). In this section, we briefly identify the anal- 
ogous symmetries described in in light of the modified non-extensive entropy S u . 
The reader is referred to [ 1 1 for background. 



(p) — q 1 ■ w 
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(i) 



9 = 1 



' Sir ^ 

, at , 



Fig. 1. Symmetry in the sequence of operations: The limit 
of the derivative equals the derivative of the limit. 



In Tsallis' original incarnation of non-extensive statistical mechanics, he used his 
"Type 1" constraint J^iPifi — U (see |7|) where the exponent q as in p\ was absent. 
This led to the following expression for the stationary probability 



(37) 



where Z' q is the obvious normalizing constant. Tsallis notes that the form of pi is essen- 
tially the same as that of p. L except that 1 — q replaces every occurrence of q — 1 and vice 
versa including in the exponents. Here, we simply note that solving the corresponding 
Lagrangian function in dl7> using the constraint 



U (38) 



leads, not surprisingly, to the analogous equation 



[1 - a(m l)Jl 



(,-l)(m-l) 

Pi(t) = L ' " v " ~' J ~ J - (39) 



which, since a = (q — l)/t, has the same reversal in sign observed in [ 1 1. 

Fleischer 1 1 1 further showed that systems with the constraint 1381 also entail a sim- 
ilar scale invariance as previously highlighted except that the corresponding energy 
transformation function is also modified with a sign change. In the context of S u and 
the constraint in (I38> . the energy transformation function that is the basis of the scale 
invariance is given by 

p m— 1 

7 m- 1 _ h (40) 

1 - a(m - l)/™" 1 



where a = (q — l)/t as before. Again, the analogous sign change in the denominator 
of J40t versus the denominator in (I25> is present. 
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5.1 Scale Invariance Using Other Constraints 

Using the same approach as in 1 1 1 and in Section it follows that for all aggregated 
energy levels A, and using the constraint in (I38> leads to the scale invariant form 



dt 



PA 

t 2 



/jr 1 - a™- 1 ) 



(41) 



where again the only difference between this and the earlier result is that every occur- 
rence of pi and f i is replaced with a pi and fi, respectively. 

Second Moments Proceeding in the same fashion as in ( 13 1 1 through fl36i and using 
analogous definitions {i.e., a 2 corresponds to the variance of the values of/™" 1 ) we 
obtain the result 

3(/ m - 1 ) _ (p m - 2 ) - a™- 1 } 2 + (i - u)(f 2m - 2 ) 



dt 



t 2 

a 2 n , (1- U )(/ 2 — 2 ) 



r- 



t 2 



(42) 



(the symbol Q serves as a reminder that these values are based on the variation over the 
entire landscape). Scale invariance in second moments with the Tsallis Type 1 constraint 
is indicated by 



dt 



t 2 

~o\ (l-u)(f 2 ™- 2 ) j 

t 2 



(43) 



where 1421 and (I43> are similar to (I34> and (I36i except that, as before, every occurrence 
of u — 1 and / has been replaced with al-ti and /, respectively. 

5.2 Probabilities and Energy Relations 

A number of additional forms of symmetry relating to the constraints are described in 
1 1 1 and their forms in the context of S u are easily inferred. Thus, e.g., defining fi in 
terms of fi we have 

fi = (44) 



1 - a(m - 1)/, 



m — 1 



where again we use a — (^r^) for notational convenience. Using this we can define 
the probability pi in terms of fi and obtain 



Pi(t) 



l-a(m- l)fi 



rn — 1 



(g-l)(m-l) 



(45) 
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which as in [T] is exactly the same form as in J37i above where the Type 1 constraint 
was used except with the /j replacing the It follows therefore that 



w , w 



Z u = £[l + a(m l)/™- 1 ] V=™*=V = £ [l - o(m - l)//"^ 



(,_i)( m _i) 



(46) 



The following lemma expands on Lemma 1 in 1 1 1 
Lemma 1. For a// u > 1. 



^K//"" 1 = "X''-/'" ' = Zl^U™- 1 )- (47) 

iVoo/: 

It follows from the definition of p,, that for all i, 

( \\ + (H-l) /•™-l]T J ^\ U 

Pif™ = ( - — L - L ^ — J — 1 /r _1 - ^ 

Now observe that j^— = — 1. Consequently for all i, 



Z^zr 1 [1 + (a=l) 

K 

7U—1 fi i / 1 \ -fin— 11 ' 



Substituting this into (I48> and simplifying we get 

p?/™ -1 = ft/r -1 (49) 

and summing over all i the result follows. ■ 

We further note that other analogous symmetries as in [ 1 1 are present with regard to 
the relationships involving the Type 1 and Type 2 constraints. These relationships sim- 
ply highlight the significance of the energy transformation function which the next sec- 
tion explores further. The following delineates these relationships the proofs of which 
correspond to those in 1 1 1 and so are omitted here. 

Lemma 2. For any energy index i and parameters a > and m > 1, define the k + 1"' 
iterate of the energy transformation of fi in terms of the k"' iterate by 

rm—l 

em— I _ {itk 

l + a(m-l)/™ 

with fi,o = fi and 

; fi 



[l + a(m-l)/ i " 
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Then, 



kk = — ■ (50) 

[1 + ka(m - I)/™" 1 ] 



Proof: 

This proof is a straight forward extension of the one in (Q. I 



5.3 A Exponential Form of Powerlaws 

In 1 1 1, Fleischer notes how the energy landscape transformation leads to an exponential 
form for a power-law distribution. It was shown that for 

<r1 

(51) 



1 + ax~< 
then 



-Xx 



C x ~ C 2 x~^ (52) 



where C\ and C 2 are constants. The exponent 7 in ( 15 11 1 served to generalize the basic 
form of the landscape transformation function so that the power law exponent could it- 
self be more general. This 7 was a component that was added for this specific purpose. 
With the definition of S u and the related energy transformation function in (125b . how- 
ever, this "artificial" addition becomes unnecessary as the exponent m — 1 in d38i more 
naturally yields the general exponent of the related power law. Consequently, the con- 
straint in (I38i and the attendant definition in d25t yield a general form for the power law 
using an exponential form based on first principles. The following theorem is therefore 
stated and the proof is omitted as it again follows analogously from the one in ffl . 

Theorem 2. Let a > and x > be such that a(m — l)a; m_1 > 1 and define 

™m— 1 



1 + a{m - l)x r ' 



using the energy transformation function defined earlier. Then for all a > 0, A > and 
m > 1 

e- X£ - d ~ C 2 x 1 - m 



as x 



00 where the constants Ci — e x / a i m 1 )andC 2 ^a 2 (m-l) 2 Ae A Mm-i)_ 



Proof: 

This proof is a straight forward extension of the one in 1 1 1. Thus, by substituting m — 1 
for every occurrence of 7 and a(m — 1) for every occurrence of a in [1J, the result 
follows. ■ 
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5.4 The Derivative of /» 

Finally, another important characteristic of the energy transformation function is worth 
noting and explains why we have used the exponent m — 1 in the previous expressions. 

Lemma 3. For any a > ( this is mathematically equivalent to the condition that q > 
1) and m > 1, 

da dq 

Proof: 

The proof is straight-forward using the chain rule. Thus, taking the derivative, we get 

d_( fi 

On 



[l + aim-l)^- 1 ]- 



[l + a(m-l)/ t 

£m 

J i 



m— 11 



[l + a(m-l)fr 1 } m - 1+ ^ 

-ST 



[l + atm-l)/™- 1 ]- 1 



zsr = ~fl 



It is apparent from the foregoing, that the energy transformation function in \25l is the 
solution of the differential equation in stated in Lemma[3]with the following boundary 
conditions: for all rn > 1, a = => fi = f. 



6 Discussion 

Fleischer 1 1 1 suggested that the scale invariance properties associated with the aggre- 
gation of energy states may provide additional perspectives on macroscopic power-law 
behavior. The appearance of the energy landscape transformation in these abstract ag- 
gregations and its relevance to power law distributions indicated in Theorems ^ and [2] 
and the more generalized entropy form S u suggest further interesting connections to 
complex systems theory which are briefly discussed below. 



6.1 Complex Systems as Graphs 

The foregoing scale invariant properties are based on aggregating energy levels. Thus, 
portions of the energy spectrum denoted by a label A are lumped together and con- 
sidered as having an energy value of The energy transformation functions also 
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suggest that as these "energy levels" aggregate they lose energy at a rate and in propor- 
tion to the magnitude of the energy level associated with the aggregation. Thus, if one 
can make a connection between aggregations of energy levels and the aggregations of 
'systems', a number of new possibilities arise in modeling complex systems. 

This notion is based on the fact that it is the components of systems that either 
contribute to or detract from the energy spectrum of a larger, more complex system. 
Thus, aggregating portions of the energy spectrum is, in some sense, equivalent to ag- 
gregating those entities that contribute to the energy spectrum. Using this perspective, 
the energy landscape transformation functions suggests that these energy-contributing 
entities themselves undergo some sort of transformation^. Thus, the energy spectrum as- 
sociated with a complex system can be modeled as a discretized energy spectrum and 
done in such a way that each such aggregation has energy values disjoint from other 
such aggregations. Depending on the level of discretization, it may even be possible 
for each such aggregation to have nearly identical energy spectrums yet not have any 
energy levels in common. Consequently, one can define subsystems of a larger complex 
system as nodes in a graph where each node possesses an energy spectrum. Arcs can 
then be used to relate these subsystems to each other in any convenient and useful way. 

Modeling complex systems using graph theoretic means is of course not new (see 
1 8 1 for a fuller account). But the scale invariant properties described above seem to fa- 
cilitate a more direct connection between statistical mechanical relationships and gen- 
eral modeling approaches for complex systems. For instance, it would be interesting 
to examine the implications of this aggregation concept in the context of Ising spin 
glass models where there are aggregations of lattice points or apply it in the context of 
Markov random fields. 

Another interesting application is in modeling systems in a far-from-equilibrium 
condition, something readily accomplished using the power law aspects associate with 
the energy transformation functions. In systems that are far-from-equilibrium, the 'local 
temperature' is not constant but "fluctuating on a relatively large time scale or spatial 
scale [9 p. 220]." But the notion of temperature for general, complex systems is some- 
what problematic. How does one define 'temperature' in systems that do not involve 
large ensembles of particles or "thermodynamic" behavior yet exhibits something akin 
to being in a far-from-equilibrium condition? See 1 9 Ch.7] for a complete discussion 
on the concept of 'temperature' in statistical mechanics. 

One recent approach for dealing with these issues has been through the use of su- 
perstatistics. Sornette describes this approach: 

A particular class of more general statistics relevant for nonequilibrium sys- 
tems, containing Tsallis statistics [] as a special case, has been termed 'super- 
statistics' [] A superstatistics arises out of the superposition of two statistics, 
namely one described by ordinary Boltzmann factors e~P E and another one 
given by the probability distribution of (3 |9 p. 220]. 

We believe this concept of 'superstatistics' is readily captured through the use of the 
foregoing energy landscape transformations and application Lemma [2] that concerns 
some recursive properties of the energy transformation functions. 

These recursive relationships, first indicated in 1 1 1 and generalized here, show that 
in terms of the energy landscape, there is a certain equivalence between values of q, 
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the number of times k a system or subsystem "undergoes" an energy transformation, 
and the temperature t. Recall, that the energy transformation comes from an expression 
involving the aggregation of energy values. But the recursive aspects of described in 
Lemma |2] suggests that a larger number of transformations and hence aggregations is, 
in some sense, equivalent to higher values of q and/or lower values of t and vice versa. 
For example, given a particular /j (or for that matter, Ja), a three-fold transformation 
fi t 3 with q = 2, m = 2 and t = 1 yields 

/ i ,3( ? = 2 > t = l) = T A_ ! 

which is equivalent to the value of a single transformation when q = 4, to = 2, and 
t = 1. When q — 2, to = 2 and t = 1/3, however, fi i also yields f+hr. ■ Thus, 
successive aggregations at some given temperature are equivalent, in some sense, to 
fewer aggregations at a lower temperature because = 1) = fi t \{t = 1/3). 



7 Conclusion 

This article proposed a slight modification of Tsallis' entropy formulation. This mod- 
ification retains the basic form of Tsallis entropy yet provides for greater flexibility in 
modeling energy constraints. The additional parameter, to, can thus be used as an ex- 
ponent of the energy value in the energy constraint equation. Doing this leads to a more 
generalized energy transformation function that also retains all of the scale invariance 
and symmetry properties highlighted in 1 1 1. Moreover, this more general energy trans- 
formation function yields a more general (and natural) power law expression given in 
exponential form. 

A number of possible issues were also examined in interpreting these results. Fu- 
ture research examines how these energy transformation functions can be leveraged 
in a number of applications, specifically, in developing new approaches for Monte 
Carlo Markov Chain simulation methodologies. Hopefully, the properties highlighted 
here will also improve our understanding of and capability to manage large, complex 
systems — something that appears to be an increasing challenge in our time. 
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